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Abstract. We analyse the low-temperature behaviour of the Heisenberg model on a 
two-dimensional lattice of finite size. Presence of a residual magnetisation in a finite- 
size system enables us to use the spin wave approximation, which is known to give 
reliable results for the XY model at low temperatures T. For the system considered, 
we find that the spin-spin correlation function decays as l/r''( T ) for large separations 
r bringing about presence of a quasi-long-range ordering. We give analytic estimates 
for the exponent rj(T) in different regimes and support our findings by Monte Carlo 
simulations of the model on lattices of different sizes at different temperatures. 
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The long history of the Heisenberg model in two dimensions is characterised by a 
competition between two contrary opinions about properties of this model. The early 
observation, made by Peierls pQ, about long- wavelength lattice waves destroying the 
localisation of particles on their lattice sites in two-dimensional crystals was followed 
later by a similar result for 2D magnets of continuous symmetry where the spontaneous 
magnetisation is destroyed by long-wavelength spin waves. Proven mathematically 
by Mermin and Wagner [2], this fact denied the very possibility of a ferromagnetic 
phase transition in this type of systems. Although the high-temperature series for the 
Heisenberg and XY models in 2D, presented by Stanley and Kaplan [3 J approximately at 
the same time, gave indication of a phase transition in both models. Being qualitatively 
similar in those works, the two models had quite different developments afterwards. 
The 2D XY model has become famous for the Berezinskii-Kosterlitz-Thouless (BKT) 
transition [1] to a quasi-long-range ordered (QLRO) phase. This special type of ordering 
cannot be characterised by an order parameter in an infinite system and manifests 
in a power law decay of the spin-spin correlation function with distance. The low- 
temperature properties and critical behaviour of the XY model on a 2D lattice are 
governed by interactions between topological defects which appear in the system [5j[6j[7]. 
This scenario is deeply connected to the symmetry of the model. In the XY model 
rotations of a spin form an Abelian group what allows formation of stable topological 
defects like spin vortices and others, in this sense it can be called an Abelian model 
in contrast to the non- Abelian ones. The Heisenberg model is non- Abelian, this is the 
main reason to deny a possibility of a BKT transition in it, since no stable topological 
defects (instantons) can be formed. The crucial evidence for an absence of a phase 
transition in the 2D Heisenberg model came from the renormalization treatment made 
by Polyakov [HI [9]. Now it is commonly believed that this model does not exhibit any 
phase transition at non-zero temperatures, although there are still some controversies 
(see e.g. Refs. [TQ], [TTl [T2J [T3], ) and reports were made that disagree with the outcome 
of Polyakov's work, assuming the possibility of a phase transition at finite temperature 
in the 2D Heisenberg model, similar to the BKT transition in the 2D XY model (see [7] 
and references therein). The last question is important in the context of an asymptotic 
freedom of QCD at AD M\ ■ 

The discussion above concerns infinite systems. And it is an infinite 2D model of 
continuous symmetry for which the Mermin- Wagner theorem has been proven. However, 
either in Monte Carlo simulations or even in reality we always deal with finite physical 
systems. It is well known now that in a finite 2D XY-spin system below the BKT 
transition temperature there is a non- vanishing magnetisation which tends to zero only 
in the thermodynamic limit [T6J [T7] . This observation goes back to Berezinskii himself 
and is supported by experimental measurements (see Ref. [18]). 

Although there is still no definite answer for the question whether the 2D Heisenberg 
model can pass to a QLRO phase or not, it is reasonable to assume that this model 
considered on a finite lattice will certainly possess some ordering, i. e. non-vanishing 
magnetisation, at low temperatures similar as the 2D XY model does. This assumption 
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is clear from the obvious fact that in a finite system, transition to the ordered phase at 
T = (when all spins of the Heisenberg model are pointed in the same direction) must 
be continuous. Hence we must see an appearence of some ordering as the temperature 
approaches zero. This is confirmed by MC simulations on the Heisenberg model in two 
dimensions [19|. 

Due to the above reasonings, we assume all spins S r = (S*, S%, S%) in the 
Hamiltonian of the Heisenberg model on a two-dimensional lattice with spacing a and 
sites defined by a vector r: 

H = ~\ E E J ( r " r ') 6 s ? 5 ? + s y r s y r , + s z r s z r ,) , (i) 

r r' 

being pointed approximately in the same direction at low enough temperatures. This 
allows us to treat the model by means of the spin- wave approximation (SWA). We 
consider the case of the nearest neighbours interaction potential J(r — r') = J5| r _ r /| )0 , 
and 1/2 stands in (PQ) to prevent double count of each bond. 




Figure 1. The angle variables 9r and 9? that can be used to define the position of 
the spin S r placed at the site r. 

We choose a special system of the angle coordinates tf'', 9^ (see Fig.l), defined 
by the relations: 



S x r = 
qy _ 

s z r = 



cos^cosfl^ , 
sin^cosfl^ , 
sin^ 2 ) , 



(2) 
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with —7i < QV) < 7i, — | < 9^ < |. Note that the variables chosen are just slightly 
modified angles (p, 9 of the spherical coordinates: 9^ = (p, 9^ = 9 — tt/2. Assuming 
that the angles 9^ , 8^ are small at low temperatures for all spins of the system, this 
choice of coordinates enables us to make use of the SWA and to substitute the scalar 
product of two spins that stands in 

S*S y r , + S y r S y r , + S z r S z r , = cos(#« - 9$ ] ) cos(0f _ 0<?>) 

+ (l - cos(0« - 0<p)\ sin 0< 2 > sin 9® , (3) 

by an expression quadratic in ^r 1 '*, 9^\ Thus, fl3]) can be written in the SWA as 

SfSI + SfS y 2 + S{S Z 2 « 1 - - 9 2 1) ) 2 - i(^ 2) - 9 { 2 2) ) 2 , (4) 

and the Hamiltonian ([T]) is reduced to 

H = H + H XY ({9^}) + H XY ({9^}) , (5) 

where 

= i E E J ( r - r ') & - (6) 



is the Hamiltonian of the 2D XY model on the same lattice taken in the SWA [20J . Hq 
can be regarded shift in the energy scale. 

The spin-spin correlation function with the assumption about smallness of all 9^\ 
9t reads: 

G 2 (R) = <S r • S r+R ) « (cos(^) - ei%) cos(0?> - 9i%)) , (7) 
where the angular brackets stand for the thermodynamic averaging: 
(...) = -Tr(...e"^), with Z = Tre"^ 

z 

and 

Tr ... = J[ i- J * d0® cos 0® .... (8) 

In the 2D XY model the power law decay of the spin-spin correlation function with 
distance serves as an indication of QLRO. A suitable quantity to characterise the decay 
of the correlation function with increase of the distance R is the temperature dependent 
exponent: 

\nG 2 (R) 

V(T) = - hm . (9) 

R^co In R 

In the case of the 2D XY model the SWA gives for the exponent n XY [20\: 

n XY = 1/(2tt/?J) , (10) 
that is reliable for small temperatures 
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Due to the separation of the angle variables 8^, 0^ in (jHJ) and ([7]) we can write 
for the correlation function: 

G 2 {R) = G^ (R) x (R) , (11) 

where 

(#\R) = ±-(4x)-» (J] £ e'^ Y ^ 008(61?) - (12) 



Gf{R) = Y 2 (U j^\d6?)cos6^ e-^« e(2) »cos(^-^ R ) .(13) 



and 

f-vr/2 
-7T/2 

and Z 2 respectively originate from the integration over 6^ and 6^ in the partition 
function Z = Z\Z 2 . Although it may be believed that the SWA applied to 0{n) models 
automatically leads to (n — l)n XY exponent, the presence of the cosine in the integration 
element in ( Fl3l) makes the problem more involved. 

Now, to define the decay of the spin-spin correlation function G 2 (R) for large 
distances R it is enough to find the asymptotic behaviour of G ( 2 1 \R) and G 2 (R) in 
the limit R/a — > oo. It is easy to see that G <y 2 \R) ) Eq.( Tl2l . is just the correlation 
function of the 2D XY model, G 2 Y (R), the asymptotic behaviour of which is well 
known: 

G { 2 1 \R)=G 2 XY (R) » (R/a)-« XY (14) 

with n XY given by ( FlOi) . So, the problem is to evaluate G 2 \R). We will follow the same 
scheme that has been used to find G 2 Y (R) [20]. For this purpose we pass to the Fourier 
variables: 

= -4 £ ^ e * = 4? E e_ikr ^ (15) 



k i" 

where iV is the total number of lattice sites and k spans the 1st Brillouin zone. Then 
the Hamiltonian ([6]) reads 

H XY ({9}) = </$>xM-k 

k^O 

with 7k = 2 — cos fc^a — cos k y a, and the integration in (fl3|) must be changed to 

^r^ k ) (i6) 




with the product taken over a half of the 1st Brillouin zone denoted by B/2. The cosine 
cos(# - 0<5r) in (USD 

can be presented in the Fourier variables as 



Re exp 



V k^O 
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with 7/£ = coskr — cosk(r + R) and r/£ = — (sinkr — sink(r + R)). The product of 
cosines in ( TT3T) can be expressed in the SWA as 



ncos^^n 



e 2 1 r ' = exp 



(17) 



using the equality ^2 r (0^) 2 = Xlk^o ^ k ^~ k we obtain an integrable form for G^\r). 
After simple integration we have 



sin 2 k 



Recall that we are interested in the long-distance behaviour of the pair correlation 
function (JTj) at low T. We estimate the asymptotic behaviour of (fT8l) in the limit 

R/a oo, /3J — > oo for two cases: 

° 2 W \ (fl/a)"^ f 0r M « 1 . (19) 

Thus, it is either constant with respect to R or equivalent to f|T^|) depending on the 
asymptotic value of the ratio ^ffi . 

Substituting (TH?j) into (ITTj) with the known expression for G^\R)-, Eq.flTJ), we get 
the following asymptotic behaviour of the spin-spin correlation function: 

G 2 {R)~{R/a)^ XY for ^! « 1 , (20) 

G a (i2) ~ (R/ar XY for ^! » 1 . (21) 

Let us clarify which of the two above estimates corresponds to the behaviour 
observable in practice. An approach used in the above derivations was based on the 
assumptions about smallness of the temperature and finiteness of the lattice. Therefore, 
the limit (3 J — > oo is physically grounded, because the lower temperature we consider, 
the more likely ordering in the system is. But the limit R/a — >• oo as well as iV — ► oo 
used to obtain approximate estimates of the integrals in practice is limited by finiteness 
of the lattice: (R/a) 2 < N. As the lattice size grows to infinity the ordering disappears 
and our approach may become invalid. From the above arguments we conclude that for 
a system of a finite size in the low-temperature limit the estimate (1201) holds. Moreover, 
power-law asymptotics (120|) . (|2~T1) brings about QLRO present in the system. Recall that 
these formulas were obtained by means of the SWA. Applicability of the latter to 2D 
Heisenberg model has been justified at the beginning of this paper by a finite system size 
that leads to residual magnetisation at low T. Therefore, although the QLRO was not 
considered as an intrinsic property of a finite-size system, the power-law asymptotics 
( !20|) . (I2~T|) gives arguments in favour of its presence. 
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Let us quote another argument that favours the QLRO in the 2D Heisenberg model, 
which by no means is connected to the SWA and holds also for an infinite system. In 
Refs. [Hj arguments were given that the model undergoes a freezing transition at non- 
zero temperature [12] with typical low-temperature configurations in a form of a gas 
of superinstantons [l1 1] . Subsequently, an onset of the low-temperature QLRO phase is 
characterised by the power-law decay of the pair correlation function with an exponent 
rj = r] XY [131 EH]- Note, that our result (120]) does not contradict the estimate of [13J, 
since the former is valid at low temperatures, whereas the latter holds at the transition 
temperature. Although our approach can not give a definite answer about the presence 
of QLRO phase in an infinite system, this question is still under discussion. 




o-o L = 8 
□—a L = 16 
A — A L = 32 
v— v L = 64 
<X> L = 128 
x— x L = 256 

Figure 2. Monte Carlo results for XY (to illustrate the well established case) and 
Heisenberg models for finite two-dimensional lattices of different sizes. 

To verify our analytic results we have performed Monte Carlo simulations of the 
Heisenberg spin model on lattices of different sizes and at different temperatures. 
The Wolff's cluster algorithm was used for this purpose [22J. The exponent rj is 
obtained on the base of three different observables, analysing the finite-size scaling of 
the magnetisation, M ~ L~2 V ( T \ the pair correlation function, G2(L/2) ~ (L)~ v ^ 
and the magnetic susceptibility, \ ~ L 2 ~ v( - T \ All three quantities are computed 
at different temperatures for varying system sizes, giving access to a temperature- 
dependent exponent f]{T). Power-law scaling found for all three quantities M, G2, 
and x supports the presence of a QLRO phase found by analytic considerations. Note 
that the lattice size in our simulations changes from A r = 8x8toA r = 256 x 256 for 
each fixed temperature (Fig. [2]) in order to obtain the finite-size scaling estimates of 
r](T) shown in Fig. [3l We cover the range of temperatures from 10~ 9 to the order of 
1. Thus in fact we observe both cases ^ R Jfl <C 1 and 3> 1. However, comparing 
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the analytic results and the outcome of our Monte Carlo simulations in Fig.3, we see 
that 1] = 2n XY fits the Monte Carlo data over the whole range of temperatures except 
the last several points (at the high temperature side of the window shown in Fig. [3]) 
which must indicate a transition to a non algebraic behaviour (see Refs. [23] and |13j). 
Thus, rj defined by 02 lj) have not been observed in our computer simulations. The 
natural conclusion is that when the condition ^g) 3> 1 is reached the temperature is 
not low enough to use the approach based on the SWA and possibly is already close 
to the conjectured transition temperature of the model (see Refs. [12] )• Hence, the 



question about a possibility to observe (!2T!) in MC simulations remains opened. But 
the important conclusion of the work is that the result rj = 2r] XY for the Heisenberg 
model in two dimensions is reliable in a wide range of low temperatures up to lattices 
256 x 256. 
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Figure 3. Comparison between the exponents r\ of the Heisenberg model obtained 
from the Monte Carlo simulations and from the analytic calculation in the SWA. The 
dashed line presents r\ XY . The inset shows an increase of the scale to make the different 
symbols used more visible. 
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